Exact Quantum-Statistical Dynamics of Time-Dependent Generalized Oscillators 
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Using linear invariant operators in a constructive way we find the most general thermal density 
operator and Wigner function for time-dependent generalized oscillators. The general Wigner func- 
tion has five free parameters and describes the thermal Wigner function about a classical trajectory 
in phase space. The contour of the Wigner function depicts an elliptical orbit with a constant area 
moving about the classical trajectory, whose eccentricity determines the squeezing of the initial 
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A quantum system of time-dependent oscillators has 
been a continuing issue of interest since the advent of 
quantum mechanics. Paul trap is one of such oscillators, 
which has a time-periodic frequency 0. Recently geo- 
metric phase has been studied for time-dependent quan- 
tum oscillators [0. Various methods have been applied 
to time-dependent quantum oscillators in many areas || . 
Agarwal and Kumar, and Aliaga et al. Q studied statis- 
tical properties of time-dependent oscillators. Also the 
density matrix and density operator for time-dependent 
oscillators were studied in Refs. || |6|. 

On the other hand, Lewis and Riesenfeld (?]] introduced 
a method to find the exact quantum states for the time- 
dependent Schrodinger equation. In particular, for time- 
dependent oscillators they found a quadratic invariant 
operator, satisfying the quantum Liouville-von Neumann 
equation, whose eigenstates provide the exact quantum 
states up to time-dependent phase factors. Even for time- 
dependent generalized oscillators each complex solution 
to the classical equation of motion leads to a pair of in- 
variant operators, linear in position and momentum op- 
erators ||. 

In this paper, using the linear invariant operators, we 
find in a constructive way the most general thermal den- 
sity operator and Wigner function up to the quadratic 
order in position and momentum operators. This den- 
sity operator is a squeezed and displaced state of a ther- 
mal one. Further, the density matrix is the thermal one 
shifted by a real classical solution, which has five free pa- 
rameters. The contour of the Wigner function follows an 
elliptical orbit with a constant area whose center moves 
and principal axes rotate along a classical trajectory. The 
shape of the ellipse measured by eccentricity determines 
the squeezing of the initial vacuum. 

The time-dependent generalized quantum oscillator is 
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described by the Hamiltonian 
H(t) 



(1) 



where X, Y and Z explicitly depend on time. Lewis and 
Riesenfeld have shown that the invariant operator satis- 
fying the quantum Liouville-von Neumann equation 



m^- t i{t) + [l(t),H(t)] = o, 



(2) 



provides the exact quantum states of the time-dependent 
Schrodinger equation as its eigenstates up to time- 
dependent phase factors. Following Ref. [|| we introduce 
a pair of linear invariant operators 



a u (t) = 
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u*{t)~Y{t)u*{t)]q 
u(t) - Y(t)u(t)]q , (3) 



where u is a complex solution to the classical equation of 
motion 




XZ-Y 1 



XY - XY 
X 



= 0. (4) 



with overdots denoting the derivative with respect to t. 
Normalizing the complex solution to satisfy the Wron- 
skian condition 



Wr{u*, u} — — (uu* — u*u) 
X 



(5) 



one can make the invariant operators satisfy the standard 
commutation relation 

[a„(t),a+(t)] = i. 

Another complex solution v to Eq. (Q), which can be 
expressed as a linear superposition of u: 

v(t) = fi*u(t) — v*u*{t), 



2 



for complex constants /i and v given by 

^ = iWr{u, v*}, v =iWr{u*,v*}, 

leads to another set of the invariant operators a v and . 
The Wronskian condition on v 

Wr{v*,v} = i & \fi\ 2 - \is\ 2 = 1, 

also guarantees the standard commutation relation 

[a v (t),al(t)] = l. 

In fact, these two sets of invariant operators are related 
through the Bogoliubov transformation 

a v (t) = p,a u (t) + val(t), 
at(t) = »*ai(t) + v*a u (t). 

The Bogoliubov transformation is written as the similar- 
ity transform ||J [Io| 

a v (t) = S^WautfSit), at(t) = §~ l (t)Si(t)§(t), 



by the squeezing operator 
S(t) 



e* tf " a >" exp 



leW'-'J cosh^ 1 - H, 



Hence this implies that by allowing all the complex it's 
satisfying both Eqs. (0) and (0) the invariant operator 
(0) is general enough for our purpose. From now on we 
shall work on the Fock bases a u and for all the complex 
u's and drop the subscript u. 

Since the invariant operator (0) satisfies Eq. (En, we 
use it to define the density operator || 



-/3X„(t) 



(8) 



Here j3 is a free parameter that may be identified with 
the inverse temperature of the system in equilibrium, and 
Z = Tr(e~^ Iu ). The density operator has five free pa- 
rameters, i.e. (3 or ujq, a complex constant 5, which is 
related with the classical position q c and momentum p c 
as will be shown below, and \fi\ and 6 V in choosing u. By 
introducing the displacement operator 

D{z) = e -*a t (*)+**«(*) j 

with z = —d/(JuJo), e = \5\ 2 / {hujo), we write the density 
operator as 



(9) 



where 



hit) 



-I3hu> a) (t)o(t) 



where /i = \/i\e ie ^,is = \is\e ". We may use the freedom 
in choosing the overall constant phase of a v , which is not 
physically important, to fix the phase 0^ = JhJ. Thus 
there are only two parameters and 9 V or a complex 
constant is, i.e. \is\ and U . 

The most general, quadratic, Hcrmitian invariant op- 
erator constructed from the pair a v and takes the form 

iv(t) = ^al 2 (t) + ^[at(t)a v (t) + a„(t)fit(t)] + ^a 2 v (t) 



+Dat(t)+D*a v {t)+E, 



(6) 



where A and D are complex constants, and B and E are 
real constants. By choosing /i and is, i.e. u such that 

A[i* 2 + 2Bfi*is + A*v 2 = 0, 

the invariant operator (^) can be written in the canonical 
form 

J u (t) = hu aUt)a a (t) + Sal(t) + S*a u (t) + e. (7) 
where 

huj = Ap*is* + B(\fi\ 2 + \is\ 2 ) + A* [iv, 
6 = Dfi* + D*is, 

e = E + -tiuio. 



is a thermal density operator. It follows that Z = Zt due 
to the unitary transformation (^). The coherent state, 
defined as a(t)\z, t) = z\z, t), is also given by 

\z,t)=D\z)\Q,t), 

where |0,t) is the vacuum state that is annihilated by 
a(t). The position and momentum expectation value 
with respect to the coherent state is 

(z,t\q\z,t) = VK(uz + u* z*) = q C) 

(z,t\p\z,t) = -^q c + ^(iiz + u*z*) =p c . (10) 

The q c and p c satisfy the classical Hamilton equations 

q c = Xp c + Yq c , 
Pc = -Yp c - Zq c . 

Now, from the definition of the thermal expectation 
value 

(6) = Tr[6p(t)] = Tv[b(z)6b\z)fc], 

we find the expectation value of position and momentum 
operators 

(q) = 1c, (P) =Pc; 
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and that of quadratic operators 

(q 2 ) = Hu*u(l + 2n) + q : 
(?) " 



(^(qp + pq)) 



where 



— (u* - Yu*){ii - Yu){l + 2n)+p 2 c , 
X z 

-^[(u* - Yu*)u + (ii - Yu)u*] 
x(l + 2n)+q c p c . 



1 



is the mean number density of Bose-Einstein distribution. 
The vacuum result is obtained by taking the limit /3 — > 
oo. It is worth noting that the dispersion of position 
and momentum around the classical trajectory (q c ,p c ) is 
entirely determined by the thermal one: ((q — q c ) 2 ) = 
(q 2 ) T and ((p-p c ) 2 ) = (p 2 )t, where (tf) T = Tr[6p T (t)]. 
Using t)(z) = e ^qa/2h eiqa f,/h e ^ Pc q/h^ we find the CQ _ 

ordinate representation of the displacement operator 

(q\D{z)\q') = e ^/ 2h e-^'/ h S(q + q c - q'), 

and that of its Hermitian conjugate, (q\D' (z^q 1 ) = 
(q'\D(z)\q)* . Hence the density matrix is given by 

P(q,q') = (qW)W) 

= (q\D\z) J d qi \ qi )( qi \p T J dq 2 \q 2 )(q 2 \b(z)\q') 

= e^-^ h p T (q-q c , q >-q c ,t), (11) 

where is the density matrix for the thermal state 
given, for instance, in Ref. ||. Finally, the Wigner func- 
tion is given by 

1 f°° 

P(q,p) = -r / dyp(q-y,q + y)e 2i ^ H 

= P T (q-q C: p-p c ), (12) 
where Pr is the Wigner function for the thermal state: 

, , tanh(/3Sw /2) 2 tanh(/3^ /2) 

Pr{q,P) = — exp 



We(<?,p) 



ujqU u \ p - 




The Wigner function (^2|) also has five parameters: q c , 
p c , f3 (or oj ), n, and 9 V . 

The Wigner functions P and Pt and their vacuum 
limit are positive definite in contrast with those for ex- 
cited states that may take negative values in some re- 
gion of phase space Hence P and Pj may be 
used as a distribution of phase space for the quantum 
evolution. For instance, the harmonic oscillator with 
X = 1/X,Y = and Z = rriQUj 2 , has the solution 



u{t) 



_ p —iu t 



/y / 2moWQ which recovers the well-known 



Wigner function with We(<Z,p) = H(q,p). In general, 
He depicts an ellipse centered at the origin, which can 
be written in the canonical form 

He = (\+p 2 + X-q 2 ) x ^, 



u u + 




u*u ( dh\(u*u) 1 / 2 \ 2 



d.t 



1 u u / 

7u + ~\ 



u*u /dhi(u*u) 



dt 



1/2 



where (q, p) are new phase-space coordinates rotated with 
the angle 



tan[20(i)] = 



2(u*u/X)[d\u(u*u) 1 ' 2 /dt] 



u*u+ (l/4u*u) - (u*ulX 2 )[d\n{u*uY/ 2 /dt}- 



Since A+A_ = 1, the area of the ellipse does not depend 
on the solution u. Therefore, as shown in Fig. 1, the con- 
tour of the Wigner function ([l2]) follows an elliptical orbit 
with a constant area whose center (q c ,Pc) in turn moves 
on a classical trajectory and principal axes (q, p) rotate 
with the angle 6(t). The shape of the ellipse somehow 
determines the squeezing of the initial vacuum |Q . 

Now we introduce a geometric measure for the squeez- 
ing of the initial vacuum and particle production in terms 
of the eccentricity of the ellipse. For that purpose we con- 
sider the exactly solvable oscillator 



tnatpi 



X 



1 



Y = 0, Z = m[uj\ - u 2 tanh(t/r)]. 



The oscillator has an asymptotic frequency u>i — {uj\ + 
Wg) 1 / 2 at t = — oo and uif = (uj 2 — w 2 ) 1 / 2 at t = oo. 
As r is an interval for the frequency change, the adi- 
abatic change is prescribed by the condition r ^> 1. 
The solution that has the correct asymptotic form u = 



I \J2muJi &t t = — oo is given by 



u(t) 



\J2mio. 



•■2Fi(-i~(u)i + uj f ), -i~(u>i 



<"/);! 



ITUJi 



where 2 F\ is the hypergeometric function. At t = — oo 
we obtain He = P 2 /(2m) + muj 2 q 2 /2 = Hi. Whereas, at 
t = oo the solution has another asymptotic form 



u(t) 



1 



where 
a±(r) 



\/2muji 

r(l - iu}jT)T(^iufT) 

r(i- i f(^± W/ ))rH§( Wl ± W/ ))- 

oo, I or | — * and u(t 



For the adiabatic limit r 
oo) = a + e~ lUft /y / 2muji. As |a+(r = oo)| = y/uJi/ojf, 
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we get He — p 2 /{2m) + mw%q 2 /2 = Hf. The contour 
of the Wigner function at t — ~ oo is given by an ellipse 
centered at the classical position and momentum 

q c = qi cos(wii + fi), Pc = -muiiqi sin(wjt + <Pi), 



where 



2h 



z , e 



and at t = oo, centered at 
q c = q f cos(w/t + iff), p c = -rmjfqf sm(u ft + <ff), 



where 



90 



2h 



\/iz + i/*z*\, e" 



lip 



fJLZ + V* z* 



mujQ \fiz + v*z*\ 

In fact, the contour is similar to the epicycle of an ellip- 
tical orbit moving on another elliptical orbit. Now the 
eccentricity e(t) = y/ min{ A± (t ) } / max{ A± (t) } and the 
rotation angle 6(t) of the small elliptical orbit explicitly 
depend on time through the magnitude of the complex 
solution 



Hi>*e- 2iuot + p*ve 2luJot . 



where 



Qf 



2h 

rruvi ' 



\a+z\ 



a+z 
\a+z\ 



From the Bogoliubov transformation between the invari- 
ant operators (||) for Hi and Hf follows the complex 
squeezing parameter 

v = im{u*(t)u}(t)-u}(t)u*(t)] 




In terms of the eccentricity e, and e/ of Hi and Hf we 
obtain the expression 

|(1 _ e?) l/2_ (1 _ e 2)l/2| 



2[(l- e f)(l- e pi/4 • 

Thus we show that the geometric shape determined by 
the eccentricity of the ellipse measures the squeezing of 
the vacuum state and thereby the amount of particle pro- 
duction. 

As an illustrative but nontrivial application of our gen- 
eral Wigner function, we consider the harmonic oscillator 
with 



X = 



1 



Hence the density matrix and Wigner function for v ^ 
describe various kinds of nontrivial quantum states be- 
yond the standard static ones. 

In summary, we showed that the most general, 
quadratic, Hermitian invariant operator (|^) can be writ- 
ten in the canonical form (Q) by suitably choosing a 
complex solution u to the classical equation of motion 
([|). Then the general thermal density operator takes 
the form (^J), which is nothing but the displaced state 
of the thermal state (||). The density matrix ( |ll|) and 
Wigner function (p"2|), the main results of this paper, are 
the thermal ones shifted by a classical solution. The gen- 
eral Wigner function for thermal state has five free pa- 
rameters, which characterize the classical configuration 
(q c ,Pc) together with (3 (or u>q) for thermal equilibrium, 
(or \v\) for the squeezing of the initial vacuum, and 
9 U . Also it is positive definite and describes a distri- 
bution of phase space for quantum evolution. Further, 
we showed that the contour of the Wigner function de- 
picts an elliptical orbit centered at a classical trajectory. 
Our general Wigner function provides a nontrivial time- 
dependent Wigner function even for a time-independent 
harmonic oscillator, whose contour describes an analog 
of epicycle of an elliptical orbit moving about another 
elliptical orbit. 



Y = 0, Z = tow, 



which has the most general solution 
1 



u(t) 



y/2mwo 



As shown in Fig. 1, the contour of the Wigner function 
moves on a small elliptical orbit about another elliptical 
orbit for the classical trajectory 

q c = q cos(u t + (/J ), Pc = -muj q sin(uj t + tp ), 
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FIG. 1: The contour of the Wigner function P. The ellipse 
centered at the origin describes the classical trajectory, a pe- 
riodic motion, and those centered on the ellipse correspond 
to a thermal state. The <p is the polar angle of the classical 
position (q c ,p c ). 



